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ABSTRACT 

We consider the reconstruction of the CMB lensing potential and its power spectrum of the full sphere in presence of sky-cuts due 
to point sources and Galactic contaminations. Those two effects are treated separately. Small regions contaminated by point sources 
are filled in using Gaussian constrained realizations. The Galactic plane is simply masked using an apodized mask before lensing 
reconstruction. This algorithm recovers the power spectrum of the lensing potential with no significant bias. 
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1. Introduction 

The gravitational potential of the large scale structure gener- 
ates a deflection of the trajectories of the Cosmic Microwave 
Background (CMB) photons. This effect, known as CMB lens- 



ing has long been predicted (Blanchard & Schneider 1987 
|Bernardeau|[T998} |Zaldarriaga & Seljak||1999| see |Lewis & 



Challinor (2006]) for a review). It was first detected by cross- 
correlating CMB data with large-scale structure ( |Smith et al. 
|2007 Hirata et al. 2008 ) and was recently detected internally us- 
ing only CMB observations ( |Das et al.|201 1| |van Engelen et al. 
|2012| l. CMB lensing consists in a remapping of the underlying 
unlensed temperature field. The lensed CMB is related to the 
unlensed by the deflection angle d, which is the gradient of the 
lensing potential (d(n) = V0(n)). The lensing potential </>(n) is 
defined as the integral along the line of sight of the gravitational 
potential of the large scale structure weighted by a geometrical 
kernel 



<p(h) 



I 



drj 



ArOteckO?) 



T(*A,i7), 



(1) 



where *P is the Newtonian potential, 77 is the conformal time, 
77 rec is the epoch of last scattering, and x is the angular diame- 
ter distance in comoving coordinates. The lensing potential thus 
contains valuable information on the large-scale structure and on 
the cosmological parameters governing their growth such as the 
equation of state of Dark Energy or the sum of neutrino masses 
( Lesgourgues etal.|2006| [Benoit-Levy et al.|2012| i. It is therefore 
of utmost importance to be able to recover the lensing potential 
and its power spectrum from the wealth of high quality CMB 
data which will soon be available, especially with the forthcom- 
ing data of the Planck satellite. 

CMB lensing modifies the Gaussian structure of the primary 
anisotropies and generates a correlation between the temperature 



and its gradient (Hu 2000). These couplings can efficiently be 
used to construct an estimator, quadratic in the observed temper- 
ature, that can be applied to data to recover the lensing potential 
dOkamoto & Hu|2003| ). 

Reconstruction of the lensing potential on the full sky re- 
quires specially designed techniques to account for the compli- 
cations inherent to real-data analysis. In particular, the presence 
of regions of the sky that are too contaminated by Galactic emis- 
sion or by point sources prevents the direct and straightforward 
use of the quadratic estimator. Various techniques were devel- 
oped to account for sky -cuts in a full sky CMB lensing analysis. 
|Perotto et al. ( 2010| ) developed an algorithm to reconstruct the 
Gaussian fluctuations of the CMB inside the masked regions. 
Smith et al. ( 2001) chose to perform the optimal filtering of the 
observed temperature, by inverting the full variance of the data. 
This expensive operation, which requires the inversion of a ma- 
trix (Npix x Npix) would be feasible at the Planck resolution but 



requires CPU consuming algorithms. Finally, |Plaszczynski et al. 
( 2012[ ) recently proposed an hybrid method that decomposes the 
region outside the Galactic mask in small patches on which a 
flat-sky analysis can be performed. 

Is this article, we present a new and simple method to re- 
construct the lensing potential on the sphere. We treat differ- 
ently the zones of the sky which must be removed because of 
the presence of point sources and the regions contaminated by 
the Galaxy. The pixels corresponding to point sources are in- 
painted with Gaussian constrained realizations. The idea is to 
restore isotropic Gaussian CMB fluctuations at the expense of 
adding a little noise in the data, so that the lensing estimator has 
the same response as it would have in the ideal case. Regarding 
the Galactic cut, we choose to apodize the Galactic mask and 
perform the usual isotropic filtering required by the lensing esti- 
mator on the masked data. As we will show later, this operation 
is almost harmless to the lensing signal. 
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In Sec. [2] we present the formalism used in this paper for 
lensing reconstruction. Then we treat each aspect separately. In 
Sec. 3 we investigate the question of the point sources mask 
and deal with the Galactic mask in Sec. 4. Finally, in Sec. 5 we 
present a reconstruction where the two effects are present. 

2. Lensing reconstruction 

We review in this section the formalism of CMB lensing and 
CMB lensing reconstruction. 

2.1. Notations 

Throughout the paper, unlensed quantities are denoted with a 
tilde. Thus, Q denote the lensed TT power spectrum and C{, the 
unlensed. CMB lensing consists in a remapping of the unlensed 
temperature on the sky: 



0[n] = 0[n + d], 



(2) 



where d is the deflection angle. This deflection angle can be ex- 
pressed as the gradient of the lensing potential 0, but for gen- 
erality can also have a curl component ( [Hirata & Se ljak 2003 1 
|Namikawaet al.|2012| i 



d = V0 + V X l/r. 



(3) 



At first-order approximation, the deflection is a pure gradient 
and the curl component is null (tfe = 0), but there exist other 
mechanisms that could lead to non-zero curl in the deflection 
angle (e.g., primordial background of gravitational waves from 
inflation ( |Cooray et al.|2()05)). Given the smallness of the signal 
that can be expected in the curl reconstruction we consider in the 
following that tfe = 0. 

2.2. Quadratic estimator 

The idea of lensing reconstruction is to rely on the off-diagonal 
terms in the CMB covariance that are created by CMB lensing. 
Lensing reconstruction then simply consists in computing the 
correlations between a filtered version of the observed tempera- 
ture with a filtered version of its gradient (Okamoto & Hu 2003). 
The first step in the reconstruction consists in filtering the ob- 
served temperature by its variance. This can be done in several 
ways, either by inverting the full variance (S + N) (Smith et al. 
2007 ), or, simply by applying an isotropic filtering. We chose the 
latter option and our inverse filtered map then takes the follow- 
ing form 



dim - 



nobs 

Im 



c! 



(4) 



where 0* s are the harmonic multipoles of the observed temper- 
ature, and C' ot the total power spectrum of the observed temper- 
ature. This total spectrum can be decomposed as 



Cf = Cf + N e , 



where C^ A is a fiducial lensed power spectrum and Ne the power 
spectrum of the instrumental noise, which is considered here as 
white and homogeneous 



N { = AT 2 e« M V. 



(6) 



Here AT is the instrumental sensitivity and 0fwvm is the angular 
resolution parameter of the CMB experiment ( |Knox|1995] l. 



The inverse variance filtered multipoles are then used as a in- 
put to the quadratic estimator. We follow the construction of the 
quadratic estimator by Okamoto & Hu (2003 1, but we slightly 
modify it along the lines of |Lewis et aL[ ( |201 Instead of filter- 
ing the inverse variance filtered map by the unlensed temperature 
power spectrum, we use the lensed power spectrum. This small 
modification has been shown to greatly simplify the expression 
of the variance of the lensing estimator ( |Hanson et al.|20TT) - 
hence the estimation of the lensing power spectrum. To summa- 
rize, we consider the following fields 

(1) (n) = Yj MM&Ytm<fi, (2) (fi) = 2 MMgYtmm, (7) 



I'm 



with 



fiiO = CeMO = 



C? 



(8) 



We t hen consid er the product of (1) with the gradient of (2) , 
(1) V0® ( |Hu|2000| . This quantity is a vector field from which we 
can extract the curl-free g( m and gradient-free C( m components. 
After calculations, ge, n and c lm take the following forms: 



and 
where 



i F(li,l 1 ,m x ,m 1 J,m)ee, m fie 1 m 1 Pte,e- 1 



(9) 



(10) 



F(t u €i,m lt tm,A,/i) = (-iy/i(A)/ 2 (4) V^2 + l)x (11) 

(12) 



(2£ x + l)(2£ 2 + 1)(2A+ 1) 



4n 



-1 1 



mi 1112 -pi 



and 



Peeiei 



(13) 



ge m and ce m are two estimators of the lensing potential and 
the curl mode. The expectation of these two estimators are 



{gem) = [A*]" 

and 
{c (m ) = 0, 



(14) 



(15) 



where A P is a scalar function that renormalizes the estimator so 



V 



(5) [A*] 



that the normalized estimator, ge„, = A,ge m , is unbiased 
(gem) = 4>e m - 

We can easily show that takes the following form 

^{({+ l) YUKAWA) 



(16) 



(17) 



(2^ + mh + 1)( -1 1 ) V^2 + \)pu A 



Q, VA(A + 1) 
C t! ^h(h + 1) 



h £1 ( 
-1 1 

-1 1 
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Operationally, this quadratic estimator can easily be imple- 
mented using the HEALPix ( |G6rski et al. 2005 1 routine 
map2alm_spin with the two components of 6> UJ V# (/; as an input. 
It is worth noting that our estimator is different from that 



of 



Okamoto & Hu 



(|2003| by a fa ctor y/t{{ + 1). This can be 



explained as Okamoto & Hu (2003) considered the divergence of 

tr\ pvtro r-t the* «-*iit"l_T-rf»*a 



the temperature-gradient product 6^' V# (2) to extract the curl-free 
part and then the lensing potential. Their estimator then involves 
an additional derivative which explain the -\J£(i + 1) factor. 

The variance of the lensing potential estimator is given by 
( [Okamoto & Hu|2003||Kesden et al.|20"03")|Hanson et al.|201 1| > 



UMg UM ,) = 5 LU 5 MM> [Cf + Nf g) + N%*>] 



(18) 



where Nf is a bias term which emerges from the un- 

lensed CMB. More precisely, the A^°' bias is the response of 
the quadratic estimator to the underlying unlensed CMB which 
has a non-zero unconnected trispectrum ( |Hu|200l) 



N 



(0,g) 



71 t,c 2 



(19) 



(2* 1 + 1)(2* 2 + 1)(* % f)c™0^ 2 

/l(A)/ 2 (^2)V^2 + l)( o % \\ + 

/i(W 2 (A)VA(A + d( o -ill- 

The N*- term can be expressed in the full-sky formalism 
(Hanson et al. 2011 1, but we need to resort to the fiat-sky ap- 
proximation dKesd en et al.||2003) > f or a computab l e expr ession 
( |Lesgourgues et al.|2005| . Following |Hanson et al.| ( |201 1| , jV (1 ' g) 
takes the following form 



N 



(i,g) 



2AI 



(2L + 1) ^f^L+X) 

fhLt 2 



(20) 



2 (-1)** 



h L 

ICfC^ 2C%C% \ A h L' 



t\ £2 

Cfffe, u h fhL'u Pie, e 2 Pie 3 tj>vt& P + L't 2 i, 



with 

fhLh 

and 

F( lL ( 2 



Cf,F< 



CiLt : 



+ Cf,F, 



y/L(L+l)i 2 (£ 2 + l)^ 



-1 1 



(22) 



{2£i + \)(2£ 2 + 1)(2L + 1) 



An 

This expression for Ft 1 u 1 i s more general than the one frequently 
used in the literature and is valid whatever is the parity of the 
quantity l\ + £ 2 + L. 

We also consider the variance of the estimator of the curl 
modes. Similarly, we have 



{clmcI m .) = 5 LU S MM , [Nf c) + N?- c) ] . 

e similar to A^ '* 
ence being in the parity conditions in their definitions 



where Nf' c) + JV[ U) are similar to Nf' g) + N ( ^ g) , the only differ- 



N 



<0,c) 



An 



(24) 



(2A + iK2^ 2+ i)^ 1 \ [Yr:cf Pu 

/l(W2«!)V^ + 1 )(o -1 l) + 

/i(W2(A)VA(A + i)(5 % {] 



N' 



(l.c) 



2AI 



(2L + 1) y[UL + 1) 

ftii& 



(25) 



2 (-l/ 2+ ^ 



ft 



2Cfc™ 2C°ie°i I u 



i\ i 2 



C L ,fi,L'h fh l> u P Ltl e 2 Pu, uPue^Pt' 



It has been assumed in Eqs. (20i and (25i that the power 
spectrum of the curl modes is null. It should be noted that the 
^y(i,c) t ernii despite being related to the curl-like part of the de- 
flection depends on the lensing potential power spectrum. This 
means than CMB lensing will be a contaminant for any future 
studies that aim at detecting faint signals in the curl modes of 
the deflection angle. Here, we follow Cooray et al. ( |2005 1 and 
use the curl modes in the lensing reconstruction mainly as a sys- 
tematic test as we expect a null signal, except for the small N^'^ 
term (see also |van Engelen et al.| ( |2012"] > for a fiat-sky discussion). 

Given the expression of the variance of the quadratic estima- 
tor, we can construct an unbiased estimator of the spectrum of 
the lensing potential ( Kesden et al.|2003) 



L £ ~ KgLMg LM ) ~ N L -N L . 



(26) 



This estimator is unbiased and at second order in <p its variance 
is given by 



(27) 



Similarly, the estimator of the power spectra of the curl modes 
reads 

cf = (c L Mcl M )-N^-N^\ (28) 

The curl term is null at second order in in the expansion of the 
remapping equation. Besides it can be shown that computations 
of the next order yield a result compatible with 0. We then have 



(Cf > = 0, 



(23) an d me variance is given by 



(29) 



(30) 
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2.3. Simulations 



In order to test the validity of the lensing estimator in the vanilla 
case of a pure lensed CMB with a perfectly white and homo- 
geneous noise on the full sky, but also in presence of sky-cuts, 
we generate simulations of lensed temperature field. We use the 
algorithm presented in Basak et al. (2008). Our lensed tempera- 
ture maps are then filtered in harmonic space using a Gaussian 
circular beam with Ofwhm = 5 arcmin to account for the instru- 
mental beam. We then add to these maps a AT = 50/vK. arcmin 
homogeneous white noise. The resulting synthetic maps have 
characteristics similar to those that the Planck experiment is ex- 
pected to produce ( |Planck HFI Core Team et al.|201 l|l. Th is map 
will be pixellized using the HEALPix ( Gorski et al.|2005| l pack- 
age. Given the experiment characteristics, we will only use the 
Nside = 1024 resolution. 

As a benchmark of both our implementation of the lensing 
quadratic estimator and the quality of our simulations, we first 
try to reconstruct the lensing signal on those perfect full sky 
synthetic sky realisation. Results are presented in Fig.[T] where 
we show the average over N = 300 reconstructions on lensed 
maps. Thorough examination of the reconstruction for the gradi- 
ent (middle panel) and the curl (bottom panel) modes reveals the 
presence of small residuals. Those are small compared to second 
order corrections (N) ' g ). They can either come from higher or- 
der corrections (that we are ignoring) or lack of convergence of 
our estimate as we are only using 300 realizations. In the follow- 
ing we will investigate how this residual degrades when taking 
into account and correct for masks. 

Now that we have demonstrated that in the ideal case of an 
uncut sky the use of the quadratic estimator leads to an unbi- 
ased reconstruction, we successively address the questions of the 
point source mask and Galactic mask. 



3. Treatment of the point sources mask 

In the following results, we use a realistic point source mask con- 
structed from the Planck Early Release Compact Source Catalog 
(Planck Collaboration et al. 2011). More precisely, we consid- 
ered the compact sources from the 100, 143 and 217 GHz chan- 
nels and masked disks with radius equal to three times the values 
of the beam of the corresponding channel. We therefore use a 
point sources mask which is composed of holes of various sizes, 
about 30, 21 and 15 arcmin. Some of the holes overlap, result- 
ing in an enlarged distribution of holes sizes (Fig. [2}. We an- 
ticipate on the following section, and we will only include in 
this mask the point sources which are outside a realistic Galactic 
mask which masks about 20% of the sky. 



3. 1 . Filling the point source mask 

The first step of the lensing reconstruction pipeline consists in 
restoring a (fake) signal in the region contaminated by resolved 
point sources. The aim of this operation is to restore the Gaussian 
statistics of the temperature map to ensure that the response of 
the quadratic estimator to this inpainted map will be unchanged 
compared to the full sky case (i.e. without masking). We will 
do so by filling the point source holes by a Gaussian realization 
constrained by the signal around the masked region. This will 
add a little bit of noise in a small area of the sky (« 2%, given 
the point source mask we use). We will see that this is a small 
price to be paid for the simplicity gained this way. 
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Fig. 1. Reconstructed power spectrum of the lensing potential 
averaged over 300 full-sky lensed simulations. Red dots repre- 
sent the gradient mode of the deflection (i.e. the lensing poten- 
tial). Green dots represent the reconstructed curl modes. Middle 
and bottom panels represent the residuals over 300 lensed simu- 
lations of the gradient (middle) and curl (bottom), with the cor- 
responding first-order bias term (black lines). Histogram errors 
bars represent the theoretical dispersion expected for one single 
reconstruction (Eqs. |27| ) and ( |30) , Error bars on the data points 
are the realization variance on 300 simulations. 



We recall here quickly the equation of the Gaussian con- 
strained realizations. The observed temperature map is decom- 
posed in two parts 



(31) 



where T\ represents the regions of the sky which are masked and 
need to be inpainted and T 2 the regions of the sky which are used 
to constrain the value in the masked regions. We also write the 
covariance matrix S of the temperature field as 



£11 £12 

£21 ^22 



The joint probability of (T; , T 2 ) is given by 
1 



P(Ti,T 2 )oc 



VdetfT) 



exp 



'2 \T 2 



T 2 



(32) 



(33) 
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mean power spectrum of the unmasked maps to the input power 
spectrum. This green line, as expected, nicely varies in the region 
defined by the black lines, which represent the expected vari- 
ance for 100 realizations in the (slightly wrong at low €) gaus- 
sian approximation for the Q. We clearly see that using only 10 
bordering pixels yield to an excess of scatter which is probably 
correlated given its apparent oscillatory behavior. However us- 
ing 15 or 20 neighboring pixels seems to be acceptable and at 
the level of precision of our simulation, the power spectrum is 
nicely recovered without any significant bias. 



Fig. 2. Histogram of the sizes of the holes in the point sources 
mask. 













10pix 


15pix 













The probability of T\ knowing the constraint T 2 , is a 
Gaussian centered on 



and with a variance 



12S 22 £ 2 l 



(34) 



(35) 



Given this formalism, the generation of a local Gaussian con- 
strained realization is straightforward (Hoffman & Ribak 1991). 

We start by generating a random realization of the tempera- 
ture T with a power spectrum Cj? d . As the variance of the masked 
regions given the constraints does not depend on the values of 
the field, both T\ and t\ have the same variance. We just need 
to shift the mean of the inpainted region to recover the expected 
mean in the masked region. Specifically, we fill the masked re- 
gions with the values 



(36) 



The computation of this quantity requires the inversion of £22- 
This operation can be prohibitively expensive if we use the full 
CMB as a constraint. However, given the statistical properties 
of the CMB temperature anisotropies, the data far away from a 
given hole will have a small constraint. This is shown in Fig. [3] 
where we represent in real space the CMB correlation function 
as well as the filter we are constructing above. As we can see 
the filter is quickly decreasing. This means that we can probably 
reduce the size of the matrix S22 by on ly computing it for the 
unmasked pixel in a narrow region around each hole and get a 
reasonable approximation of the result. 

We perform tests on the Gaussian constrained realizations 
on our point sources mask with three different choices for the 
size of the constraint region around the hole, using a 10, 15 or 
20 pixel border (resp, 34.4, 51.5 and 68.7 arcminutes with our 
HEALPix map resolution). The results of these tests are shown 
in Fig. [4] which presents the ratio between the mean power spec- 
trum of 100 maps whose point source mask as been filled using 
the procedure described above and the power spectrum used to 
perform the random realizations. The different results (depend- 
ing of the size of the border region used to build the constraint) 
are to be compared with the green line showing the ratio of the 









— positive part 

— ■ negative part 




ft' 


ifi 





Fig. 3. Top panel is the CMB correlation as a function of the 
angular separation. Bottom, the filter corresponding to EuS^l 



as defined in eq. (36 1. Vertical lines show the angular separation 
corresponding to 10, 15 and 20 pixels (i.e. 34.4, 51.5 and 68.7 
arcminutes) 



3.2. Point sources and lensing reconstruction 

We now wish to demonstrate that our treatment of the point 
sources mask is harmless for lensing, i.e. that 1) it does not create 
any spurious lensing signal, and 2) that it enables an excellent re- 
construction of the power spectrum of the lensing potential. We 
first test the case of unlensed maps. In that case, the terms in- 
volving the lensing power spectrum disappear and according to 
Eqs. 26 and 28 the estimated spectrum should be equal to zero. 



Here again using around 15 to 20 bordering pixels is required 
in order to get an unbiased reconstruction. Figure [5] shows the 
mean over 100 lensing reconstructions on unlensed maps. Both 
the cases with 15 and 20 bordering pixels prove to be in very 
good agreement with the unmasked case, for both the gradient 
and curl modes. We thus confirm that the point source inpainting 
does not create any artificial signal that could be seen as related 
to lensing via the lensing estimator. 

We now turn our analysis on lensed maps. In that case, we 
need to modify the estimator of the lensing power spectrum to 
correct for the fact that the regions that are inpainted are filled 
with a perfectly Gaussian realizations and should therefore not 
contain any lensing signal. The outputs of the lensing estima- 
tor must therefore by rescaled by a factor that corresponds to the 
fraction of the signal that is masked. This factor should in princi- 
ple be scale dependent as the large scale of the lensing potential 
are hardly affected by masking sub-degree scale pixels. However 
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Fig. 4. Ratio of the mean of the N= 100 power spectra of the in- 
painted maps to the fiducial power spectrum for three values of 
the width of the region used for constraining the Gaussian real- 
ization, 10 (blue), 15 (red), and 20 (magenta) pixels. The green 
line represents the case where no mask is applied. Black lines 
are the cosmic variance expected given the number of simula- 
tions and binning. 
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Figure J6] represents the residuals of the reconstruction, i.e. 

ACt^ - ct^' PS - cf^' fld averaged over 300 simulations. We also 
show the residuals of the curl estimator. We compare those re- 
sults to the full-sky case, with no inpainting. As can be seen, 
there is no significant difference between the two cases, which 
validates the use of our inpainting algorithm for lensing recon- 
struction. 

We therefore have a tool that can be efficiently used to in- 
paint the point sources and that can recover both the 2-point 
and 4-point statistics of a Gaussian field without introducing any 
spurious signal or bias. 



H 



{ | border=15 
f } No mask 



t | border-15 
f f No mask 
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Fig. 6. Residuals of the reconstructed power spectrum of the 
lensing potential averaged over 300 full-sky lensed simula- 
tions in presence of the point sources mask. Top panel presents 
the gradient mode of the deflection (i.e. the lensing potential). 
Bottom panel presents the reconstructed curl modes. In both 
cases, we compare with the full-sky case (red for the gradient 
modes and green for the curl modes) Histogram errors bars rep- 
resent the theoretical dispersion expected for one single recon- 
struction. 



Fig. 5. Reconstructed gradient (top) and curl (bottom) modes of 
the deflection angle on 100 unlensed maps. Point sources are in- 
painted using 15 (red) to 20 (magenta) bordering pixels to con- 
strain the Gaussian realizations. The unmasked case (green) is 
shown for comparison. 



we found that using a scalar scale-independent renormalization 
gives an excellent reconstruction of the lensing power spectrum. 
The modified estimators thus reads 

c H>,ps _ (8lm8* lm ) -N[' g) _ ^ (hg) 
fps 

where / P s correspond to the fraction of the sky which is un- 
masked by point sources. With our fiducial point source mask, 
fps = 0.98. 



4. Treatment of the Galactic mask 

The second problem frequently encountered in full-sky lensing 
reconstruction is the presence of large-scale foregrounds, mainly 
along the Galactic plane, that strongly contaminates the lens- 
ing signal. As mentioned in the introduction, several possibili- 
ties exist to deal with this complication. We choose a new ap- 
proach which consists in masking the contaminated regions of 
the Galactic plane and then directly applying the lensing esti- 
mator on those masks maps. We first begin by presenting some 
analytical computations on the effect of the mask on the recon- 
structed lensing potential, and we show how these effects can be 
efficiently alleviated to provide a robust and rapid way of recon- 
structing the lensing potential. 
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4.1. Analytical computation of the mask 

In presence of a mask, the expectancy of the lensing estimator 
become 



?LM> - / l 



M 



Luff + 



IMF 



(38) 



where M LM is a coupling matrix which replaces the scalar nor- 
malization A* and 0^ is a lensing independent quantity which 
we will loosely refer to as the "mask mean field". 

The variance of the estimator applied on a masked tempera- 
ture map becomes an intricate expression 



LM T L,M 



^ L.M ^ 



MF 



(39) 



where C^ F denotes the power spectrum of the mask mean-field. 



Niff and NY^' are terms similar to the full-sky equivalent N\ 

and N^' g) but also depends on the structure of the mask. All the 
terms in the two previous equations involve high order integrals 
of functions which can be written as intricate summations of 
Wigner 3-j symbols. They are directly related to the structure 
of the mask and depend on its harmonic coefficients wg m . Those 
functions and integrals can be written down analytically, but can- 
not lead to a form which can be efficiently evaluated numerically. 
However, we find that under certain circumstances, notably de- 
pending on the mask (and we will provide precise example in 
Sec. |4.2) >, those coupling matrices are essentially diagonal, in the 
sense that their effect on power spectra is just to apply a constant 
normalization factor. This normalization factor is related to the 
fraction of the sky which is unmasked (Hivon et al.|2002| ). More 
precisely, we have 



(().,?) 



and 

N (0 ' g) * 

L,M 

where 



M'Lm 



/gal,4JV L , 



N 



(Ug) 



ft 



gal,4 



(40) 



(41) 



(42) 



The presence of the coefficients of the mask at the power of four 
is related to the fact that the variance of the lensing estimator is 
a resummation of the 4-point correlation function on the CMB. 
The variance of the lensing estimator then becomes 



(gLMg*L>M'} ~ ^LL'(>MM'fgal,4 \c\ 



] + c 



MF 
L ' 



(43) 



Even though all these mask kernels cannot be computed, 
some of them are fairly easy to estimate by a Monte-Carlo pro- 
cedure. Applying the estimator on unlensed maps which have 
the same spectral content as the lensed map will precisely give 
the mask mean field 



(8Tm) 



IMF 
"LM' 



(44) 



It is thus straightforward to construct an unbiased estimator of 
$lm by averaging over se 1 
on unlensed masked maps 



^ by averaging over several outputs of the estimator applied 



IMF 

<Plm 



1 

^MF 
uru 



(45) 



<p M¥ is an unbiased estimate of the mask mean-field and its 
variance can be approximated by 



/ ^MF IMF* \ x X 

\<Plm<Pl>m>) = °LL>0A 



£,MF,mask + /gal.4 ^ f ((),g) 



N; 



(46) 



We can then define a mask mean-field debiased estimator by 
subtracting the estimated mask mean-field to the lensing poten- 
tial reconstructed from the masked temperature field. We then 
define 



ftMF _ „ 2.MF 

8lm ~ 8lm - <P LM - 

The mean of this estimator is 



& MF 



>=z 



M 



(47) 



(48) 



and its variance becomes 



/"MF"MF*\ ^ x x t 
\SlmSl j M'' ~ °LL'OMM'ha\A 



C7 + 



l + 



(0,8) 



,(49) 



The removal of the mask mean-field by Monte-Carlo induces a 
small increase in the variance of the estimator, but this increase 
can be as small as required by averaging over a large number of 
unlensed realizations, which is a numerically cheap operation as 
it does not require the production of lensed simulations. 

The estimator of the lensing power spectrum now takes the 
following form 



/gal,4 



1 + 



1 



(50) 



It should also be noted that on the curl estimator side, the 
situation is little simpler as there is no mask-mean field (Fig. [9] 
bottom panel). Therefore there is no need to subtract a "curl" 
mask mean field to the curl estimate. 

If we now consider the variance of the estimator of the lens- 
ing power spectrum, it will depend on 



/gal,8 



— Iv 

^Vpi X 



(51) 



as it involves the 8-point correlation function (Kesden et al. 
[20031 |Hanson et al.|201 \\ . We will then have 



cr 2 (cf' mask )= fi 



J gal,4 



+ i 



M + iV (o,«) +iV (U) 



(52) 



We have set up the formalism to perform a reconstruction on 
a masked field. We will now give precise examples and explain 
the regime of validity of our hypothesis. 

4.2. Apodization 

If we apply the quadratic estimator directly to the masked tem- 
perature fields, both the mask mean field and the band couplings 
caused by the presence of the mask will highly bias the re- 
construction (Fig. [?}, both for the gradient and curl estimator. 
This pathologic behavior is caused by the sharp transition at the 
boundary of the mask, that creates important power leakage at 
all scales. 

In order to reduce the off-diagonal terms in the coupling ma- 
trix and the frequency extent of the mask mean field, we simply 
need to first apodize the mask before masking the temperature 
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— cf* ,th J $ Curl modes 

I { Gradient modes 




Fig. 7. Reconstructed lensing potential (Eq. 50 1 over 100 lensed 
simulations using the binary mask (no apodization). Green dots 
are the gradient mode which should overlap the black line, and 
magenta dots the curl modes, which should be compatible with 
zero. 



field. Apodization smoothes the spectral content of the mask and 
reduces the mask-mean field amplitude and removes almost all 
the off-diagonal terms. This operation has nevertheless a draw- 
back as it decreases the quantity of data used for reconstruction 
and therefore degrades the statistical significance of the recon- 
struction. The apodization of the mask is performed by applying 
a cos-like function to the pixels bordering the mask so that the 
apodized region goes smoothly from to 1 over the apodization 
length apo . 

Figure [8] presents the variance of the estimator of the mask 
mean field (or equivalently the estimated power spectrum of the 
mask mean-field) for both the gradient (top panel) and curl (bot- 
tom panel) estimators and for different values of the apodiza- 
tion length. In all cases, the variance has been rescaled by the 
corresponding / ga i,4 factor. In this figure, iVjS" = 100. We here 
confirm that the mask does not create a mean-field in the curl 
estimator (bottom panel). We also see that even a small apodiza- 
tion over 1 degree is enough to remove all the off-diagonal terms 
and restore the ability of the quadratic estimator to successfully 
reconstruct the lensing power spectrum. On the top panel, we 
clealy see that the mean-field frequency extent decreases when 
the apodization length is increased. We therefore choose the 
higher value, 8 apo = 10° to perform the reconstruction. With this 
value of apodization length, we have / ga i.4 = 0.65. 

The final reconstruction is presented on Fig. [9] where we 
show the mean over 100 residuals of the reconstructed lensing 
potential power spectrum. In this case, the reconstruction seems 
less in agreement with the fiducial spectrum at low multipoles, 
but it should be emphasized that the small bias that can be seen is 
well with the expected errors bars, whereas the mask-mean field 
is actually several orders of magnitude higher that the lensing 
potential at those scales. 

4.3. Normalization 

The rescaling of the power spectrum of the reconstructed map 
by 1 // g ai,4 can be thought as a modification of the normalization 
factor applied to the unnormalized estimator g( m . Indeed, rescal- 
ing the power spectrum by l// ga i.4 is equivalent to applying a 



f-l/2 



We can elaborate on this 



normalization factor equal to A f f^J A 
by estimating the normalization factor to apply to the unnormal- 
ized reconstructed potential by computing it by a Monte-Carlo 
procedure. There are two ways of doing so, either by correlat- 
ing the lensing potential reconstructed from a lensed simulation 




Fig. 8. Power spectrum of the mask mean field rescaled by / ga i 
(top, gradient modes; bottom curl modes). The full red line is the 
Gaussian noise N gc divided by the number on unlensed simu- 
lation, N=100. When the binary Galactic mask is applied (blue 
lines) the reconstruction is strongly biased at high multipoles. 
Note the absence of a mask mean-field in the curl modes. 



+ 

£ -0.4 



< 

£ 0.2 



10 30 50 100 200 300 500 

I 



Fig. 9. Residuals over 100 lensed simulations of the gradient 
(top) and curl (bottom) modes of the deflection angle in pres- 
ence of the Galactic mask, apodized over # apo = 10°. Histograms 
represent the error bars for one reconstruction (Eq. 52 1. 



with the input lensing potential used to generate the lensed tem- 
perature map, or by correlating two lensing reconstructions from 
two different temperature maps lensed with the same lensing po- 
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tential. We thus form the following correlations that define the 
Monte-Carlo normalizations Af MC ' 2 and Af MC ' 4 



1 N 



C 



(f>(f>,th 



fmVft,,) ,0,MC,2 
; = : A { 

N N 



IE 



c 



<p<p,th 



N(N - 1 ) 4^ ^ >0 tm lm (Af MQ4 ) 2 ' 



(53) 



(54) 



where , 



are to be understood as reconstructions from 



; 8tm and SL 

different realizations of temperature maps lensed with the same 
lensing potential. As Eq. ( |53"] l brings in the two-point correlation 
function of the masked temperature and Eq. ( 54 1 the four-point 



correlation function, relations between A ( 
depend on f g& \ t 2 and / ga i,4, where 



/gal,2 



1 

N~ 



AZ' MCA and Al 



(55) 



More precisely, we have 



>,MC,2 



— — — , andA^ 

/gal,2 



MCA 



(56) 



Figure [TU] presents the ratio of the normalization computed 
by Monte-Carlo to the theoretical normalization. Red dots corre- 
spond to Af MC ' 2 and green dots correspond to A^' MC ' 4 . The first 
bin is not in agreement because of the presence of the mask mean 



field which yield a non-zero correlation in Eq. ( 54 1. Accordingly, 
when computing the normalization by correlating the recon- 
structed lensing potential to the input potential, the mask mean 
field disappears from the correlation. This figure shows a very 
good agreement between the ratio and the expected constant val- 
ues, which validates once again the use of an apodized Galactic 
mask for CMB lensing reconstruction. It should be noted that the 
two normalizations computed here do not have the same utiliza- 
tion. We have seen previously that A^' MC ' 4 is used to normalize 
the power spectrum of the reconstructed lensing potential. As 
A^' MC ' 2 only brings in a single reconstructed potential in its def- 
inition, it should be used to normalize correlations between the 
reconstructed lensing potential and any other field that would 
have a non-zero correlation with the lending potential but not 
with the reconstruction noise. A perfect example for such a field 
would be an external dataset of tracers of the large scale struc- 
ture. 



n l.i- 



f,.L2 



Fig. 10. Ratio of the normalization computed by Monte-Carlo to 
the theoretical normalization in the case of a reconstruction with 
an apodized Galactic mask. 



therefore by rescaled by the Galactic mask renormalization value 
/ g ai,4- However the inpainted regions lack the lensing signal, but 
possess the full Gaussian structure, so that the Gaussian noise 
coming from those regions is the theoretical N) . The estimator 
of the lensing potential power spectrum then takes the following 
form 



^,mask+PS 



l 

■'PS 



igLMg* LM ) 



al,4 



1 + 



1 



N MF 
unl 



N' 



(O.g) 



-N 



and its variance becomes 



/gal,8 



(57) 



(58) 



PS-'gal^ 
fid 



c 



1 + 



N' 



■ N. 



a,g) 



The residuals over 300 lensed simulations of the gradient and 
curl power spectra reconstruction are presented in Fig. [TTJ We 
notice the presence of biases at low multipoles both for gradi- 
ent and curl modes and at high multipoles for the curl modes, 
but their amplitude is small, about 0.1 cr. However, the residu- 
als show excellent agreement with zero in the multipole range 
100 < I < 1000, for both the gradient and the curl modes 
of the deflection angle. This multipole range is also where the 
signal-to-noise ratio is the highest, thus validating the described 
pipeline. 



5. Point source and galactic mask 



6. Discussion 



All the ingredients are now in place to reconstruct the lensing 
potential and its power spectrum when both the point sources 
and Galactic masks are included. The algorithm is fairly sim- 
ple. We first inpaint the regions masked by the point sources 
mask using the algorithm described in Sec. 3.1 We then apply 



the apodized Galactic mask on this inpainted map. We chose 
to use the largest apodization length we considered previously, 
#apo = 10°. We then generate a set of unlensed simulations un- 
der the lensed temperature power spectrum on which we run the 
quadratic estimator to estimate the mask mean-field. 

Renormalization of the reconstructed lensing potential 
power spectrum requires some care as some of the masked re- 
gions in the Galactic plane do not have any signal and must 



We have presented a new, simple, and robust pipeline for CMB 
lensing reconstruction. This pipeline is designed to treat two im- 
portant issues in a full-sky lensing reconstruction analysis: the 
presence of point sources that need to be masked and Galactic 
foregrounds. We treat the first problem by filling in the holes of 
the temperature map by Gaussian constrained realizations. This 
operation perfectly restores the Gaussian structure of the origi- 
nal map and inpainted maps can then safely be ingested in the 
lensing estimator. 

The issue of Galactic contamination is treated by simply ap- 
plying an apodized Galactic mask to the temperature map before 
the lensing estimation. The presence of the mask creates a bias 
in the reconstructed potential. However, since this bias is solely 
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Fig. 11. Residuals over 300 simulations of the gradient (top) and 
curl (bottom) modes of the deflection angle in presence of the 
point sources mask and of the Galactic mask, apodized over 
#apo = 10°. Histograms represent the error bars for one recon- 
struction (Eq. 59 1. 



related to the mask and does not depend on the actual lensing po- 
tential, it can be efficiently estimated from unlensed simulations 
and then subtracted to the final result. This simple operation is 
enough to guarantee an unbiased reconstruction of the lensing 
power spectrum. The contribution of unresolved point sources 
to the lensed CMB trispectrum is assumed to be small and is 
therefore ignored. 

By using an isotropic filtering during the reconstruction pro- 
cess, this pipeline remains analytical, in the sense that both the 
normalization and the Gaussian noise TV*' keep their theo- 
retical expression as they are just properly rescaled by scale- 
independent factor to account for the missing power due to the 
presence of masked region. 

In this article, we have used a simple noise model by con- 
sidering that the detector noise is white and homogeneous. This 
will not be the case for a realistic experiment like Planck, where 
the scanning strategy lead to a non-uniform coverage of the sky. 
Analytical predictions of the effect of inhomogeneous noise on 
lensing reconstruction have been investigated in Hanson et al. 
(2009). The pipeline described in this work can easily be gener- 
alized to inhomogeneous noise, by taking into account the noise 
structure in the unlensed simulations when the mask mean-field 
is computed. In that case, the estimated mean-field would be the 
sum of the noise and mask mean-fields. Similarly, any other sys- 
tematic effects that would create a fake lensing potential can be 
treated in a similar way, the difficulty being more in the ability to 
correctly simulate these effects than in the lensing reconstruction 
in itself. 
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